Disorder and interaction induced pairing in the addition spectra of quantum dots 
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We have investigated numerically the electron addition spectra in quantum dots containing a 
small number (A'^ < 10) of interacting electrons, in presence of strong disorder. For a short-range 
Coulomb repulsion, we find regimes in which two successive electrons enter the dot at very close 
values of the chemical potential. In the strongly correlated regime these close additions or pairing 
are associated with electrons tunneling into distinct electron puddles within the dot. We discuss 
the tunneling rates at pairing and we argue that our results are related to a phenomenon known as 
hunching, recently observed experimentally. 

PACS numbers; 73.20.Dx, 73.40. Gk, 71.30+h, 73.20.Jc 
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In a small metallic island weakly coupled to the en- 
vironment the number of electrons is quantized at low 
temperatures. Because of the Coulomb repulsion from 
the electrons already on the island, it takes a finite en- 
ergy to add one more electron to the island. This can 
be achieved with the help of an an external gate voltage, 
coupled capacitatively to the island. Additions of single 
electrons occur roughly periodically as a function of the 
gate voltage. This is the essence of Coulomb blockade 
(CB), one the most robust facts in mesoscopic physics. 

By means of an experimental technique known as 
single-electron capacitance spectroscopy (SECS) [0,1), 
one can study electron additions from a metallic electrode 
to a semiconductor quantum dot. In these experiments 
the electrons tunnel into localized states of the dot one by 
one, starting from the very first electron; each tunneling 
event is recorded as a peak in the differential capacitance, 
measured as a function of the gate voltage. 

A few years ago a SECS experiment [yH] showed that 
electrons sometimes entered the dot in pairs rather than 
individually, thus violating all the common wisdom that 
we have on CB. A more recent and systematic investiga- 
tion [Q has shown that quite generally, in dots contain- 
ing N < 200 electrons, electron additions are not evenly 
spaced in gate voltage. Rather, in some cases, they group 
in bunches of up to 6 electrons. Experimentally the first 
bunching occurs already when the number of electrons 
N 7 and with increasing N it evolves from occurring 
randomly to periodically at about every 5th electron. To 
explain this puzzling phenomenon two different theoreti- 
cal models have been proposed |^,^ , in which two paired 
electrons enter the dot coherently. Classical simulations 
10-^ are able to reproduce some experimental features 
observed in the presence of a strong magnetic field. How- 
ever the mechanism responsible for the pairing effect is 
still unclear. 

The purpose of this work is to carry out a full quantum 
mechanical calculation of the addition spectra of dots 
containing a small number of particles < 10, including 
their spin. Our goal is to investigate the regime where 
aperiodic electron pairing first appears. We consider the 



limit of strong disorder generating localized states that 
can be occupied by two electrons of opposite spin [p^ . 
When the Coulomb repulsion is short range, we find two 
different regimes where successive additions can occur 
sequentially but anomalously close to one other. The 
first case takes place for intermediate values of the direct 
Coulomb interaction but strong on-site repulsion, which 
favors the appearance of a polarized dense droplet with 
no holes or doubly occupied states. Both electrons par- 
ticipating in the pairing tunnel into the edges of the dot 
but at spatially distinct regions. The second situation 
occurs in the strongly correlated regime, with strong val- 
ues of the direct Coulomb interaction competing with the 
on-site repulsion and the disorder. In this case pairing 
is characterized by the formation of distinct puddles of 
electrons; doubly occupied states appear at the center of 
the dot, where one the two electrons tunnels. Specific 
features of pairing in this regime allows us to relate these 
results to the bunching phenomenon of Ref . |Q] . 

We model the quantum dot as a finite portion of a 
square lattice with J\fx x JVx sites, described by the tight- 
binding Hamiltonian 

H = J2{e,-eVM,.c,,,, 

+*E (cLs,-exp'^-+H-C.) + i?int (1) 



where the e.^'s are site energies randomly distributed 
with uniform probability between —W and W; Vg is a 
gate voltage. The hopping term, connecting nearest- 
neighbors sites, is proportional to the matrix element t 
modified by Peierls phases Lpij , due to the possible pres- 
ence of a magnetic field, a —] , | is a spin variable. The 
interaction part is given by 



XI ^^^^ 

i>j ;cr,(T' 
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where ni^^ is the number operator, v{Ti — Tj) is the ma- 
trix element of the direct part of the Coulomb interaction 
and U is the on-site-interaction constant. 
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FIG. 1. Number of electrons A*' on a 3 x 4 dot vs. gate 
voltage for one particular disorder realization {W — 7.5t). 
and nearest-neighbor Coulomb interaction (V = 2t). The 
on-site-interaction strength is U = 8t (dashed line) and 
U = oo (solid line) respectively. The arrows point to the pair- 
ing between the additions 6 7 and 7 8. Inset: A'' vs Vg 
for four disorder realizations of polarized electrons {U = co), 
displaying pairing between the N = 7 and N = 8 states. 



Experimentally the Coulomb interaction is screened by 
the external gate. For simplicity we will model the in- 
teraction as a nearest-neighbor repulsion, v{ri — Vj) = 

2 

VSrj^Ti+s — ^^Vj.ri+s, where a is the lattice constant 
pT| . The important parameter describing the strength 
of the average Coulomb interaction relative to the Fermi 
energy is = l/y^Tmas, where n is the electronic den- 
sity and ub is the Bohr radius. For our lattice model 
= ^ (V/At), where v = {N)/[{Af, - l){Afy - 1)], and 

(N) is the average number of electrons on the dot. 

By using Lanczos techniques, we have diagonalized the 
Hamiltonian Eq. ([|) on a 3 x 4 dot, within the Hilbert 
subspaces corresponding to a fixed total number of elec- 
trons between 1 and 10. The ground-state energy and 
wave function, the total number of electrons N = (N) 
and site occupation (rii^^), can be obtained for different 
values of the parameters entering Eq. is controlled 

by the voltage V^. At a specific value V^~^^ given by 
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the number of electrons on the dot jumps from N to 
+ 1. Here E^ is the ground-state energy for A^ parti- 
cles at Vg = 0; /i(A^) is the chemical potential, which is 
the quantity measured experimentally 

In Fig. we plot the dot occupation A^ as a function 
Vg for one random configuration in the regime of strong 
disorder {W = 7.5t) and intermediate values of Coulomb 
interaction V — 2t, corresponding to ~ 0.5. For onsite 
repulsion U = 8t (dashed line), the two additions 6 =4« 7 
and 7 8 occur at close gate voltages; in contrast, all 
other electrons enter the dot at well-spaced voltages. 
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FIG. 2. Ground-state site occupation < Ui > for different 
values of A'^, for the dot of Fig. hi (solid line). The area of 
the circles is proportional < Ui >. Pairing occurs between 
the N — 7 and N = 8 states. The grey circles indicate the 
sites where the largest portion of the incoming electron, par- 
ticipating in the pairing, is distributed. 

Moreover, the tendency of electrons 7 and 8 to pairing is 
strongly enhanced by increasing U. This is clearly shown 
by the solid line of Fig. |^ where U — oo, corresponding 
to spin-polarized electrons. In this case the two pairing 
particles still enter the dot sequentially but almost at the 
same gate voltage. Since this result refers to one particu- 
lar disorder realization, one might think that it represent 
an extremely rare episode. Remarkably however, this is 
not the case, as the inset of Fig. |l| shows: out of 8 ran- 
dom realization that we have tried, half of them displays 
paring between the additions 6 7 and 6 =^ 8, when the 
system is polarized, i.e. U = oo. 

In order to understand how pairing takes place, we look 
at how the site occupation < rii > changes when the two 
participating particles tunnel into the dot. Fig. |^ displays 
< Tii > for the 3x4 sites of the dot, corresponding to the 
polarized case of Fig. |^. Because [/ = oo, no double occu- 
pancy is allowed and the electrons tend to form one rather 
compact puddle (see the state iV = 6 in Fig. ||). This is 
reminiscent of the effect of the exchange interaction stud- 
ied within the Hartree-Fock approximation for a dot 
in a strong magnetic field: exchange generates a local 
attraction between the electrons, causing the formation 
of a dense droplet. The 7th electron, the first involved 
in the pair, tunnels mainly into two almost empty sites 
at the right edge of this compact puddle, as shown by 
the grey circles of the N=7 state. The 8th electron, the 
second involved in the pairing, fills up two already par- 
tially filled sites on the bottom edge of the dot. Thus, 
the pairing of Fig. |l| is associated with electron additions 
into spatially distinct regions of a compact electron pud- 
dle. Because the Coulomb interaction is short range, the 
energy costs of these two additions can be almost equal. 
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FIG. 3. Number of electrons A*' on a 3 x 4 dot vs. gate 
voltage for one particular disorder realization {W = 8t) and 
nearest-neighbor Coulomb interaction (V = 9t). The two 
lines correspond to an on-site-interaction strength U — 30t 
(dashed line) and U = 29.06t (solid line) respectively. The 
arrows points to the additions 7 8 and 8 9, where pair- 
ing occurs. The two lines differ perceptibly only at pairing. 
For the system relative to the solid line, the inset displays 
fj,{N) vs the magnetic flux </> for N = 7,8, 9. 

The results shovifn above hold for intermediate values 
of the Coulomb interaction « 2t <^ < 1), much 
smaller than the strength of disorder and on-site repul- 
sion. Experimentally however, pairing occurs at larger 
values Coulomb repulsion [vg > 2). If we increase the 
strength of the direct Coulomb interaction relative to 
the on-site repulsion and the disorder, we reach the limit 
where the dot finds more advantageous to generate lo- 
cal singlets of doubly occupied localized states. Instead 
of grouping into one compact puddle, the electrons now 
can form distint puddles. Electron pairing can also take 
place in this strongly correlated regime. 

An example is shown in Fig. ^, where we plot N vs Vg 
for a strength of the direct Coulomb repulsion V = 9t 
Vs =2), close to the disorder strength, W = 8t. The 
two curves corresponds to two values of the on-site repul- 
sion, U — 30t and U = 29.06t. We focus on the additions 
7 8 and 8 9, where pairing takes place Note 
that the two curves differ only in this region, indicating 
that the on-site repulsion U plays a crucial role. We can 
again understand how this comes about by looking at the 
site occupation (rii) for the states involved in the pairing. 
This is shown in Fig. | for U = 30t. The 8th electron, 
the first involved in the pair, enters the dot at the top left 
corner of the dot. The 9th electron jumps mainly into 
one of the central sites, which becomes partially doubly 
occupied. The dot shows a tendency to create two sepa- 
rate puddles of electrons, one composed of the 5 occupied 
sites along the left edge, and the second made up by the 
3 sites of the right edge plus the doubly occupied site. 
Each electron involved in the pair tunnels into one of 
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FIG. 4. Ground-state site occupation {rii) of the 3x4 dot 
of Fig. ^ (dashed line, U = 30t). Pairing occurs between the 
states N = S and A*' = 9. Most of the second electron partic- 
ipating in the pair {N = 9) tunnels into a central site, which 
becomes partially doubly occupied, (n) = 1.70. 

these spatially distinct regions. Finally, when the 10th 
electron enters the dot, the gap between the two regions 
is filled and the dot is occupied more uniformly. This 
example supports the suggestion, borne out of the exper- 
iment, that pairing is associated with electron localiza- 
tion in distinct puddles of the dot [Q . The merging of the 
two puddles upon increasing N , corresponds to a sort of 
localization-delocalization transition. 

If U is slightly smaller, e.g. U = 28t, the roles of the 
8th and the 9th electrons are interchanged: the first elec- 
tron of the pairing tunnels into the central site, creating 
the spin-singlet; next, the second electron occupies the 
top edge site. By tuning U between 30t and 28t, we have 
found that the intermediate value U = 29.06t gives rise 
to the closest additions 7 =4' 8, 8 9. These two events 
still take place sequentially, (see Fig. ||). For this value 
of U, an analysis of (n^) reveals that, in each tunneling 
event participating in the pairing, half electron goes into 
the central site and the other half goes into the top-left 
corner site. We have checked that other disorder realiza- 
tions display a similar behavior. 

The inset of Fig. |^ displays the magnetic flux de- 
pendence of the chemical potential, Eq. ^, for the system 
corresponding to the "best pairing" {U = 29.06i). For 
such strong values of the interaction and the disorder, 
the traces of the addition spectrum are only very weakly 
dependent on (f). This implies that the pairing states 
N = 8 and iV = 9 remain close to each other for the en- 
tire range of (j). Experimentally, in bunchings occurring 
at low N < 10, the traces also stick together for the entire 
range of (f), but at the same time they rise nearly linearly, 
showing that the field dependence of the localized states 
approaches the one of the lowest Landau level. 

So far we have discussed the occurrence of pairing only 
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FIG. 5. The addition spectral function for the two parti- 
cles forming a pair as in Fig. The dashed and solid lines 
refer to the first and second particle of the pair respectively. 
The peaks indicated by the arrows are proportional to the 
tunneling rates for the A'' — > A'^ + 1 process, (a) and (b) refer 
to dashed and solid line of Fig. |^ respectively. 

from the point of view of energy balance. We conclude 
our analysis by discussing the rates at which the two 
paired electrons tunnel into the dot in the strongly cor- 
related regime > 2. To this purpose we have computed 
the addition spectral function, defined by 



Af+l 



E4. 



(4) 



where |<&^^^) are the eigenstates of the {N + l)-system 
corresponding to energy E^^^. The operator ^ 



in Eq. ^ implies that electrons tunnel coherently into all 
the sites of the dot with equal probability. The term in 
Eq. ^ involving ($^"'"^| ^ c|^|$^) is proportional to 
the rate of tunneling from the ground state of the N- to 
the [N + l)-particle system. 

Fig. 1^ displays the spectral function for the addition 
of the two pairing electrons, corresponding to the situ- 
ation of Fig. ^. The peaks pointed by the arrows are 
proportional to the tunneling rates. In Fig. |^(a) we can 
see that for U = 30t the tunneling of the first particle 
forming the pair is much smaller than the second. As 
discussed above, in this case the first particle is mainly 
localized into one corner site of the dot (see Fig. ^). The 
presence of fast tunneling in one state of the pair is also 
seen experimentally, albeit starting at strong magnetic 
field and at higher N . At low N the experiment shows 
that the tunneling rates of both particles are large for all 



magnetic fields. As shown in Fig. ||(b) the two tunneling 
rates are both large when U takes on the value that en- 
sures the closest additions. This can be understood if we 
remember that in this case the electronic densities of the 
pairing particles are distributed among several sites of 
the dot. Thus, even if the system is strongly correlated, 
the tunneling rates for the pairing states are not negligi- 
ble, in agreement with the experiment. In contrast, the 
two models in Refs. ^ and Q predict a dramatic sup- 
pressions of the tunneling rate at pairing. 

In conclusion, our numerical simulations show evidence 
of pairing in the addition spectra of disordered quantum 
dots with strong on-site repulsion. Specific features of the 
pairing observed in the strongly correlated regime > 2 
suggest that our results may be related to the analogous 
phenomenon seen experimentally in the regime of a small 
number of electrons. 
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